Quantum complex networks 
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In recent years, new algorithms and cryptographic protocols based on the laws of quantum physics 
have been designed to outperform classical communication and computation. We show that the 
quantum world also opens up new perspectives in the field of complex networks. Already the 
simplest model of a classical random network changes dramatically when extended to the quantum 
case, as we obtain a completely distinct behavior of the critical probabilities at which different 
subgraphs appear. In particular, in a network of N nodes, any quantum subgraph can be generated 
by local operations and classical communication if the entanglement between pairs of nodes scales 
as iV -2 . 



On the one hand, complex networks describe a wide 
variety of systems in nature and society, as chemical reac- 
tions in a cell, the spreading of diseases in populations or 
communications using the Internet pQ. Their study has 
traditionally been the territory of graph theory, which 
initially focused on regular graphs, and was extended to 
random graphs by the mathematicians Paul Erdos and 
Alfred Renyi in a series of seminal papers [2 , 3 , 4 in the 
1950s and 1960s. With the improvement of computing 
power and the emergence of large databases, these theo- 
retical models have become increasingly important, and 
in the past few years new properties which seem universal 
in real networks have been described, as a small-world [5] 
or a scale-free [6] behavior. 

On the other hand, quantum networks are expected 
to be developed in a near future in order to achieve, for 
instance, perfectly secure communications [7J|8]. These 
networks are based on the laws of quantum physics 
and will offer us new opportunities and phenomena as 
compared to their classical counterpart. Recently it has 
been shown that quantum phase transitions may occur 
in the entanglement properties of quantum networks 
defined on regular lattices, and that the use of joint 
strategies may be beneficial, for example, for quantum 
teleportation between nodes |9j [10]. In this work we 
introduce a simple model of complex quantum networks, 
a new class of systems that exhibit some totally un- 
expected properties. In fact we obtain a completely 
different classification of their behavior as compared to 
what one would expect from their classical counterpart. 

A classical network is mathematically represented by 
a graph, which is a pair of sets G = (V,E) where V is 
a set of N nodes (or vertices) and E is a set of L edges 
(or links) connecting two nodes. The theory of random 
graphs, aiming to tackle networks with a complex topol- 
ogy, considers graphs in which each pair of nodes i and 
j are joined by a link with probability Pij. The simplest 
and most studied model is the one where this probabil- 
ity is independent of the nodes, with p iy j = p, and the 
resulting graph is denoted Gjv,p- The construction of 
these graphs can be considered as an evolution process: 
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FIG. 1: Evolution process of a classical random graph with 
N = 10 nodes: starting from isolated nodes, we randomly 
add edges with increasing probability p, to eventually get the 
complete graph Kio for p = 1. 



starting from N isolated nodes, random edges are suc- 
cessively added and the obtained graphs correspond to 
larger and larger connection probability, see Fig. [I] One 
of the main goals of random-graph theory is to determine 
at which probability p a specific property P of a graph 
Gn, p mostly arises, as N tends to infinity. Many proper- 
ties of interest appear suddenly, i.e. there exists a critical 
probability p c (N) such that almost every graph has the 
property P if p > p c (N) and fails to have it otherwise; 
such a graph is said to by typical. For instance, the criti- 
cal probability for the appearance of a given subgraph F 
of n nodes and I edges in a typical random graph is [11] : 

Pc {N)=cN- n '\ (1) 

with c independent of N. It is instructive to look at 
the appearance of subgraphs assuming that p(N) scales 
as N z , with z G (— oo,0] a tunable parameter: as z 
increases, more and more complex subgraphs emerge, 
see Tab. [T] 

We consider now the natural extension of the previous 
scenario to a quantum context. For each pair of nodes 
we replace the probability pij by a quantum state pij of 
two qubits, one at each node. Hence every node possesses 
N — 1 qubits which are pairwise entangled with the qubits 
of the other nodes, as depicted in Fig. [2] The goal is 
then to establish maximally entangled states within cer- 
tain subsets of nodes by using protocols that consist of lo- 
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TABLE I: Some critical probabilities at which a given sub- 
graph F appears in classical random graphs of N nodes con- 
nected with probability p ~ N z : cycles and trees of all orders 
appear at z = — 1, while complete subgraphs (of order four or 
more) appear at a higher connection probability [T]. 



be 




FIG. 2: An example of a quantum random graph on five 
nodes. (A) Each node is in possession of four qubits which 
are entangled with qubits belonging to the other nodes. All 
the connections are identical and pure but non-maximally en- 
tangled pairs. (B) Imperfect pairs can be converted into max- 
imally entangled ones with some probability of success p (here 
p — 0.25); this strategy mimics the behavior of classical ran- 
dom graphs. 



cal operations and classical communication (LOCC) [12], 
where each node can apply different measurements on its 
quantum system and then communicates its results to the 
rest. This is in fact a very natural scenario since maxi- 
mally entangled states are the resource for most quantum 
information tasks, and thus the goal of any strategy. As 
in the classical random graphs we consider that pairs of 
qubits are identically connected, with pij = p. Further- 
more, we restrict ourselves to the simplest case where the 
states are pure, i.e. p = \<p)(ip\ , since it already leads to 
some very intriguing phenomena. We take these states 
to be: 



\<p) = y/l-p/2\00) 



(2) 



where < p < 1 measures the degree of entanglement 
of the links. As before p scales with N and we write 
\Gn,p) the corresponding quantum random graph. The 
choice of the coefficients in Eq. Q becomes clear if 
one considers the following simple strategy where links 
are treated independently: each pair of qubits tries to 
convert its connection by LOCC into the maximally 
entangled state |$ + ) ex |00) + |11). The probability of a 
successful conversion is p, which is optimal [131 [14] , and 
therefore the task of determining the type of maximally 
entangled states remaining after these conversions can 
be exactly mapped to the classical problem. In that 
case we obtain the results of Tab. [IJ and for example 
for z > — 2 the probability to find a pair of nodes that 
share a maximally entangled state is one, whereas that 
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FIG. 3: Graphical representation of \K C ). Each node has 
exactly c — 1 orthogonal states in which c — 2 qubits are |0) 
and one is |1); qubits that were originally entangled share 
the same state. A link denotes, here, the separable state |11) 
while 1 00) is not drawn, and \K C ) is defined as the coherent 
superposition of all perfect matchings of the complete graph 
K c . In this example an appropriate labeling of the states 
leads to \K 4 ) oc |1111) + |2222) + |3333). 



of having three nodes sharing three maximally entangled 
states is zero unless z > — 1. 

Allowing strategies which entangle the qubits within 
the nodes offers new possibilities and brings powerful 
results. This is indeed a general fact in quantum 
information theory as, for example, in the context of 
distillation of entanglement [15] or in the non-additivity 
of the classical and quantum capacity of quantum 
channels [16 . As a first illustration of the advantage 
of joint actions on qubits we show how some relevant 
multipartite entangled states, like the three-dimensional 
Greenberger-Horne-Zeilinger (GHZ) [17] state of four 
particles, can be obtained when z = — 2. At each node 
we apply a generalized (or incomplete) measurement [18] 
whose elements P m are projectors onto the subspaces 
consisting of exactly m qubits in the state |1), i.e. we 
count the number m of "links" |11) attached to the 
nodes without revealing their precise location. Remark 
that such links are separable and so do not represent 
the sought connections |3> + ). For each measurement 
we get a random outcome m whose value is either or 
1 since, as in the classical case, the probability to get 
m > 2 is zero for infinite N [19 , see App. A. Setting 
p = 2cN~ 2 one finds that there are in average c nodes 
having exactly one qubit in the state |1), whereas all 
nodes where m = factor out since they are completely 
uncorrelated with the rest of the system. In Fig. [3] we 
describe the remaining state, which we write \K C ), and 
show that for c = 4 it is a GHZ state of four qutrits. 

We now turn to the main result of this paper: for N 
tending to infinity and with p ~ TV -2 one is able to obtain 
with unit probability a quantum state with the structure 
of any finite subgraph. That is, for an arbitrary subgraph 
F = ( V, E) composed of n vertices and I edges, one can 
get the state \F) consisting of I maximally entangled 
pairs |^> + ) shared among n nodes, according to E: 



\F) = ($\<S>+) Et . 



(3) 



This result implies that in the quantum case the struc- 
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ture of Tab. [T] completely changes, as all subgraphs 
already appear at z = —2. Actually, our results are 
more general since any state of some interest in quantum 
information theory, as W states, Dicke states, graph 
states, etc., also arise at z = —2. We sketch here the 
four steps of our proof and refer the reader to App. B 
for detailed calculations. First, we create the state 
\K C ) described in Fig. [5J with c = n + D, D = d 2 and 
d = 2 l : n nodes will be kept to build the final state 
while D additional nodes are needed to establish the 
desired quantum correlations. Second, we measure all 
links corresponding to the subgraph K n and the other D 
vertices are measured in the Fourier basis, which leads 
to a highly entangled state of n nodes of dimension D. 
Third, in order to extract the right correlations, we split 
each node into two subsystems of dimension <i, measure 
one of them in the Fourier basis and post-select a specific 
outcome pattern. Finally, we show that the resulting 
state, a n-partite GHZ state of dimension can be 
converted into \F) by some appropriate projections. 

At this point let us briefly describe a possible setup 
for an implementation of our ideas, where atoms store 
the quantum information, and thus represent qubits, 
while photons are used to create remote entangle- 
ment. This scenario is currently well admitted to be 
the most promising one for a realization of quantum 
networks [7]. In particular, we consider the case where 
(continuous-variable) entanglement contained in a 
two-mode squeezed light is transformed into (discrete) 
entanglement between atoms trapped in distant high- 
quality cavities [20]. Assuming perfect operations one 
can drive the system so that its steady state is exactly 
described by Eq. [2j with p < 1 due to a finite squeezing 
of the light. In App. C we discuss a simple model 
of mixed-state networks where the source of squeezed 
states fails to emit, with some probability, any light. We 
show that the quantum phenomena considered so far 
persist despite these imperfections, and we believe that 
this will still be the case for more general types of errors. 

In conclusion, we have introduced a model of com- 
plex quantum networks based on the theory of random 
graphs and have shown that allowing joint actions on 
the nodes dramatically changes one of their main prop- 
erties, namely the appearance of subgraphs according to 
the connection probability. In fact, all classical expo- 
nents collapse onto the value z = — 2 in the quantum 
case and we expect a large variety of new phenomena 
in, for instance, quantum models of small-world or scale- 
free networks. The model we have introduced is quite 
simple in the sense that connections are represented by 
pure states, whereas more realistic setups involve mixed 
states. We discussed this problem for a specific type of 
imperfections, and even if it may be much more difficult 
to tackle it in full generality we are confident of the per- 
sistence of unexpected and intriguing phenomena in real 
quantum networks. We hope that our results will inspire 



work in this direction and shed some light on the very 
active domain of complex networks. 
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Appendix A: General considerations and notation 

As in the classical case we consider the exponent z in 
p = N z to be larger than or equal to —2, since the overlap 
of the quantum random graph and the product state |0) 
of all qubits in |0) approaches unity for z < — 2 and in 
the limit of infinite lattices: 



N(N-l) 



K^|0>| 2 =(l-|) 2 ^exp(- 



-N z + 2 



and therefore no local quantum operation is able to create 
entanglement between nodes in this case. Furthermore, 
the outcome distribution of the measurements P m fol- 
lows the classical degree distribution of a random graph 
Gn, p ', with p' = p/2. In fact the only information we 
get from \Gn jP ) by applying P m is the number of links 
1 11) attached to a given node, and each such link gets 
a weight p/2 in the outcome probabilities. This distri- 
bution has been shown to be well approximated by the 
one that we get if the nodes are considered to be inde- 
pendent, so that the expected number x m of outcomes 
m reads, setting p = 2cN z and for z < — 1: 



E(x m ) 



ml (N — 1 — m)\ 



N-l- 



_j\jm(z+l) + l 



ml 



(4) 



In the text we often mention the act of measuring a node, 
or more generally of a system of dimension d: if no ba- 
sis is specified we mean a measurement in the compu- 
tational basis, i.e. we project the system onto the states 
1 1 ),..., \d). We also use the expression "to measure a 
link" to indicate that one of its qubits is measured in the 
{ |0), |1)} basis. Finally, we introduce the notation 

1 d 

|^) = ^=^e^|jj...j), fcG{l,...,d}, (5) 
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for the Fourier transform of = The latter 

state is referred in the text as a GHZ state of dimension d 
(on n nodes), and a measurement in the Fourier basis of 
a system of dimension d is its projection onto the states 

\^ d ). 



Appendix B: Proof of the main result 

We give here the details of the construction of a state 
\F), starting from \Gn,p) an d using LOCC only. But 
let us first note that if one is able to find a construction 
which succeeds with a strictly positive probability, say 
Pf, that does not depend on N but on the subgraph 
F only, then \F) can be obtained with a probability 
arbitrary close to one. The reason is that we can always 
subdivide the TV nodes into L sets of N/L nodes, with 
L ^> 1/pf, and apply the same construction on each set. 
These sets can be treated as independent if we initially 
discard all links connecting different sets, which is done 
by measuring the corresponding qubits. Note that in 
what follows we do not try to optimize the procedure 
since in the limit of infinite N we are mainly interested 
in the existence of a finite probability pp, not in its 
maximal possible value. 

First step. We start the construction by creating the 
state \K C ), with c = n + D, D = d 2 and d = 2\ which 
can be obtained with a probability approaching unity in 
the limit of infinite N. In fact one can choose a value 
d ^> c, apply the projections P m on the nodes of the 
quantum network and get a system of average size d . 
The number of nodes of the resulting state can then be 
decreased in a deterministic fashion: one measures all 
qubits of a node, gets for all outcomes except for a 
random one whose neighbor automatically factors out. 
Hence, the system is projected onto \K c /-2) and the 
procedure can be iterated until we get \K C ). 

Second step. We remove all connections shared be- 
tween n nodes of \K C ), i.e. we measure the n(n — l)/2 
concerned links, and the operation is successful if all out- 
comes are 0. In that way we build a state that is the 
coherent superposition, as described in Fig. [3j of all per- 
fect matchings of the join graph of and of the empty 
graph on n nodes. We further measure these D nodes, 
but this time in the Fourier basis {|$^' c_1 ),l < k < 
c — 1} in order not to reveal where the links |11) lie. 
We can correct the possibly introduced phases and the 
resulting state reads: 

D 

\Vn ) Yl 1^2... in)- (6) 

Third step. It is not convenient to deal with sums 
whose indices are subject to constraints, so we develop 
Eq.[6]to let the sums freely run from 1 to D. For example, 
the state on three nodes is expressed as: 

D D D 

lvf)= {\iij)+\iji)+\jii))+2^\iii). 

i,j,k=l hj = l i=l 

More generally this leads to a weighted and symmetric 
superposition of states of the form ®[ =1 \$\.) for all 
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FIG. 4: Scheme for entanglement generation between two 
nodes: two distant cavities are simultaneously driven by a 
common source S of squeezed light, and in the steady state 
of the system the two atoms A and B become entangled. 



partitions (Ai, A2, . . . , A r ) of n, where Y^i=i A* = n (a 
partition of a positive integer is a way of writing it as a 
sum of positive integers). We want to remove all terms 
of this sum but the last one, which is a GHZ state that 
will allow us to obtain \F). To that purpose we use the 
fact that = |$f ) = \§ d n ) m for all n and d, split 

each node into two subsystems of dimension d and mea- 
sure one of them in the Fourier basis. This operation is 
successful if all outcomes are k = 1, except the last one 
which should be k = d — n + 1. In fact, to see what hap- 
pens to a state |$^) we note that ($^|<I>4 d > <x l&^-f) 
if m > 1 and Sj^ otherwise. Therefore, by sequentially 
measuring (^>{' d | , a state |3>^) shared among any m < n 
nodes transforms as: 

lO = \$ d m d ) - \* d ~±i d ) - • • ■ - \^- m+14 ) - 0, 

as long as d > m. But this is always the case because, 
without loss of generality, the subgraph can be con- 
sidered to be connected, so that I > n — 1 and thus 
d > 2 l > 2 n ~ 1 > n > m. Hence all terms vanish except 
the GHZ state \® d n ). 

Fourth step. The last step consists in transforming 
\<$> d n ) into \F). To that end first expand Eq. [3J i.e. 
write explicitly all its terms in the computational ba- 
sis, and group the qubits according to the connections 
E. This leads to a sum of product states of the form 
\(fi,i) • • • \(fi,n) with i = 1,...,2^. Since we have cho- 
sen d = 2\ we can now apply the measurement element 

J2i=i I^m'X^I on eacn n °d e 3 °f l^nX which achieves the 
desired transformation and concludes the proof. 

Appendix C: A mixed-state scenario 

We consider here the setup described in [20] , where two 
high-quality cavities are simultaneously driven by a com- 
mon source of two-mode squeezed light, see Fig. |4j The 
continuous-variable entanglement contained in the light 
is transformed into discrete entanglement between atoms, 
and the steady state of the system is described by the 
non-maximally entangled pure state we use throughout 
the text. Let us now introduce some errors into the net- 
works, considering an imperfect source of light which fails 
to emit squeezed states with some probability e. Equiv- 
alently, this imperfect source produces the vacuum state 
1 00) with probability e so that the connections of the 



5 



quantum network are: 

p=(l-e)\<pM +e 100X001, (7) 

with \(f) defined in Eq. [2j For this error we can show 
that it is still possible to construct, in the regime z = — 2, 
mixed states which are close to arbitrary quantum sub- 
graphs. To that end let us go through the four steps 
of the construction described in the previous paragraph. 
The first step consists in creating the state \K C ) by ap- 
plying the measurements P m on the nodes of the net- 
work. These measurements are not affected by the im- 
perfect source since no extra links |11) are added: only 
the number of outcomes m = 1 slightly decreases from 
c to (1 — e)c, but one can choose a larger constant d in 
p = 2c'N~ 2 in order to get with certainty c outcomes 
Pi. For c = 4 for instance, the remaining nodes are in a 
mixture of the desired state \K A ) and some completely 
separable states: 

1 - x 3 

p K4 = x \K A )(K A \ H — ^2 I^X^I , (8) 

6 i=i 

with x = (1 — s) 2 for infinite N. Despite the presence of 
separable states, pka is useful for quantum information 
tasks since it is distillable for all e < 1, i.e. the coeffi- 
cient x can be brought arbitrarily close to unity if one 
possesses a large number of copies of it. However, in the 



regime z — — 2 it is impossible to get several copies of 
Pka on the same four nodes. But this is not a problem 
since, alternatively, we can repeat the construction 1/x 
times so that any use of \K C ) is still achieved with high 
probability. More generally, the state pK C on c nodes is 
a mixture of \K C ) and of some partially separable states, 
with x equal to (1 — e) c ^ 2 . Note that all terms appear- 
ing in pKc are also present in \K C ), but they now are 
probabilistically weighted. This structure is maintained 
throughout the construction of the quantum subgraphs 
(steps two to four of the proof), so that with probability 
x we create \F), and with probability 1 — x we get a sum 
of useless quantum states. Therefore, for all e < 1, with 
a strictly positive probability we can achieve quantum 
communication or computation that is impossible within 
a classical framework. 

In conclusion, we have discussed a possible implemen- 
tation of our model of quantum networks and shown that 
the proposed construction of quantum subgraphs is ro- 
bust against some specific imperfections, namely a noise 
|00)(00| from the light sources. Remark that this is not 
the case for errors involving terms like |01)(01| , |10)(10| 
or |11)(11| in the connections. However, in that case, 
we are confident that other quantum strategies (based 
on purification methods for instance) will still lead to in- 
triguing and powerful phenomena, thus stimulating fur- 
ther work on complex quantum networks. 
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